Bistable diode action in left-handed periodic structures 
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We study nonlinear transmission of an asymmetric multilayer structure created by alternating 
slabs of two materials with positive and negative refractive index. We demonstrate that such a 
structure exhibits passive spatially nonreciprocal transmission of electromagnetic waves, the ana- 
logue of the electronic diode. We study the properties of this left-handed diode and confirm its 
highly nonreciprocal and bistable transmittance by employing direct simulations. 



An all-optical diode is a spatially nonreciprocal device 
that allows unidirectional propagation of a signal at a 
given wavelength. In the ideal case, the diode trans- 
mission is 100% in the forward propagation, whereas 
it is much smaller or vanishes for backward (opposite) 
propagation, yielding a unitary contrast. The possibility 
of achieving highly nonreciprocal transmission through 
a passive, nonlinear device is not only novel and in- 
teresting concept, but also useful for many applica- 
tions such as optical isolation and all-optical process- 
ing. Such unidirectional action was demonstrated for sev- 
eral asymmetric structures and different nonlinear mate- 

nais jiaiainia. 

In this paper, we discuss a novel type of spatially non- 
reciprocal device based on an asymmetric multilayered 
structure created by alternating slabs of two materials, 
conventional and left-handed (LH), the latter is a mate- 
rial with both negative electric permittivity and negative 
magnetic permeability which results in an effective neg- 
ative refractive index Such multilayered structures 
are a special class of structures consiting of a sequence 
of flat lenses that provide optical cancellation of the con- 
ventional, or right-handed (RH), layers leading to either 
enhanced or suppressed transmission 0, 10, 111. 

We employ a general idea to obtain diode action by 
making the structure both asymmetric and nonlinear. 
Also, we consider a periodic superlattice, which is ex- 
pected to possess resonant properties, and thus enhanc- 
ing nonlinear effects. We consider an asymmetric super- 
lattice consisting of three stacks with four LH/RH dou- 
ble layers each [see Figs.nja,b)]. We assume that in the 
nonlinear regime, a dielectric defect layer with Kerr-type 
nonlinear response is inserted between stacks 1 and 2, 
as shown in Fig. ^b). We study wave transmission in 
both linear and nonlinear regimes by using the transfer- 
matrix method (TMM) 0] > and also by employing direct 
numerical simulations based on the pseudospectral time- 
domain (PSTD) method H3- The transfer matrix of the 
system analytically relates the electromagnetic fields on 
either side of an element 01 ■ Complex structures can be 
build up out of simple elements through multiplication of 
the respective matrices. Thin nonlinear elements can be 
modelled, e.g., as delta- functions prLIT^|. 

In the PSTD method, the Maxwell equations are dis- 
cretized in time and space. The spatial derivatives are 



Stack 1 



Stack 2 



Stack 3 



(a) 



(b) 



2 3 
Position [X 



o J 



FIG. 1: Schematic of the LH/RH superlattice for transmis- 
sion in (a) linear regime and (b) nonlinear diode regime. Lay- 
ers of LH material are shown in light gray. The dielectric 
defect layer is shown in dark gray. 



approximated using discrete Fourier transforms and the 
temporal derivatives using central differences. From 
this, one derives update equations for the different fields 
and by iteration propagates the fields through time. 
The problem of the inherent periodicity of the Fourier 
transform can be removed through the use of perfectly- 
matched-layer absorbing boundary-conditions pH , Il5| . 
The PSTD method samples all field components at the 
same location in the unit cell, which is advantageous 
when both the electric permittivity and the magnetic per- 
meability vary with position |lo| . An instantaneous Kerr 
nonlinear material is included by directly solving the cu- 
bic equation for the intensity |17| at each iteration. 

We model the LH material as having Lorentzian fre- 
quency dependence in both electric permittivity e r and 
magnetic permeability /i r , 



s r (u) = 1 + 



Lir(u) = 1 + 



pe 



pm 



(1) 

(2) 



Here, u> 



pe j ^pm 



are the corresponding plasma frequencies, 
wi e ,a->i m are the resonance frequencies, and 7 e ,7m are 
absorption parameters. We use u> pe = 1.1543 x 10 11 s _1 , 
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FIG. 2: Transmission spectrum of the structure and its com- 
ponents in the linear regime, (a) Stack of eight periods with 
individual layer thickness Ao/3 (solid) [Stack 1 and Stack 3 
in Fig. I^a)] and a stack of four periods with individual layer 
thickness Ao/5 (dashed) [Stack 2 in Fig. |T[a)]. (b) Combi- 
nation of these stacks when arranged as in Fig. [Ra,). The 
shading indicates the location of the band gaps for Stack 1 
(light grey) and Stack 2 (dark grey) in the limit of an infinite 
number of periods. 



uj pm = 1.6324 x 10 11 s _1 , oJu = uj lm = 2ir x 5 x 10 6 s _1 , 
and include small losses, j e = 2 x 7 m = 2ir x 6 x 10 6 s _1 . 
With these parameters, the refractive index of the LH 
material n f» — 1 at frequency /o = 15 GHz, and the 
material is left-handed for frequencies / < 18.5 GHz and 
right-handed for / > 26 GHz. The RH layer is assumed 
to be air. 

The structure shown in Fig. (Ha) consists of three 
stacks with four LH/RH double layers each. The outer 
stacks (Stack 1 and Stack 3) have individual layer thick- 
nesses of Ao/5, while the central stack (Stack 2) has in- 
dividual layer thickness Ao/3; Ao being the free-space 
wavelength of incoming radiation of the frequency /o- 
Figure ^b) shows the same system with a structural de- 
fect introduced between Stack 1 and Stack 2. The defect 
layer has thickness 2Ao/25. 

First, we study the transmission properties of the in- 
dividual components of our superlattice in order to an- 
alyze the effect of the layer thickness on the positions 
of the band gaps and the transmission spectrum. Fig- 
ure Ufa) shows the transmission spectrum of a structure 
made of the combination of Stack 1 and Stack 3 directly 
attached to each other (solid), as well as, the transmis- 
sion spectrum of Stack 2 (dashed). The spectra show 
several transmission bands and gaps, and we clearly see 
a shift of the conventional band gaps due to the scal- 
ing of the structure. However, the band gap associated 
with vanishing average refractive index (n) = (around 
fo = 15 GHz) remains at the same frequency and only 
changes its width [To|. 

The transmission spectrum of the combined structure 
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FIG. 3: Transmission spectrum of the structure in Fig. ffl b) 
for incidence from the (a) left and (b) right, for different out- 
put amplitudes, calculated by TMM. The arrows indicate the 
direction in which the peaks shift with increased input (solid) 
and output (dashed) intensity. 



of Fig. da) is shown in Fig. EJb). In the frequency in- 
tervals where the central stack is opaque, the structure 
has a band gap, while at frequencies where the outer 
stacks are opaque but the central stack is transparent the 
structure behaves similar to a cavity resonator and shows 
the characteristic Fabry-Perot cavity peaks. In this sys- 
tem cavity-like behavior occurs in the frequency ranges 
of 13.5-14.3 GHz, 16.3-17.5 GHz, and 25.8-29.9 GHz. 

To study the nonlinear transmission of our multi-stack 
structure, we introduce a defect layer with a linear elec- 
tric permittivity e r = 4, see Fig. [fl». For the TMM 
calculations, we place a delta- function Kerr- type non- 
linearity at the left interface of the defect layer 01 
e r [l + S(x — xo)], xo = 1.6Ao, while in the PSTD cal- 
culations the entire defect layer has a nonlinear electric 
permittivity e r (t) = 4+x*- 3 -* |E(i)| 2 , with the nonlinearity 
coefficient x^ 3 ) = +4. We study the difference in trans- 
mission when the light is incident from opposite sides of 
the structure. In the linear case, the reciprocity princi- 
ple holds as expected and the transmission coefficient is 
independent of direction of incidence, while a nonlinear 
asymmetric structure can show directional dependence 
of the transmission indicating the corresponding diode 
action effect Q, • 

First, we employ the TMM and study the steady-state 
results. In this situation it is convenient to fix the trans- 
mitted (output) field and calculate backwards to find the 
incident and reflected fields. In Fig. 0we present a spec- 
tral interval including the first two peaks (peaks A, B) 
above the (n) = band gap of Stack 2 for several fixed 
transmitted field strengths and different directions of in- 
cidence. The inclusion of the defect shifts all peaks in 
Fig-I21b) to slightly lower frequency in the linear regime. 
Due to the nonlinear properties of the defect, the fre- 
quencies of the transmission resonances vary with the 
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FIG. 4: TMM results: Bistability curves for / = 16.19 GHz 
(a) and 25.17 GHz (b), at equal relative detuning below a 
linear transmission peaks A and C. The curves show incidence 
from the left (solid) and the right (dashed). Peak C shows 
bistability for greater but not for smaller detuning than shown 
in (b) . Peak A shows bistability even for much smaller relative 
detuning. 



local electric field intensity which in turn varies with the 
output field. Since the behavior of peak B is very sim- 
ilar for the two directions of incidence, the output field 
strengths in the figure are chosen to highlight the behav- 
ior of peak A. The frequency shift of peak A occurs at 
very different output field strengths depending on the di- 
rection of incidence. As the input field is increased, all 
the transmission peaks shift to lower frequencies (solid 
arrows in Fig. The curves are presented with fixed 
output field and one notices that all peaks shift to lower 
frequencies with increased output field (dashed arrows in 
Fig. |3J) , except for peak A in the case of incidence from 
the left [Fig. 01a)], which moves to higher frequencies. 
This opposite shift in frequency means that while gener- 
ally the output field at the transmission peak increases 
with increased input field, for this peak the output field 
actually decreases as the input field is increased. Re- 
ferring to Fig. |2Ib) one notices that peak A is located 
slightly above a band gap of the central stack (Stack 2). 
With increased input intensity the nonlinear defect shifts 
the peak to lower frequency such that it moves deeper 
into the band gap and becomes increasingly suppressed. 
Behavior similar to peak A is also found for the first peak 
(peak C) above the first Bragg band gap of Stack 2 (not 
shown) . 

Next, in Fig. we compare the output intensity as 
a function of input intensity for incidence from the two 
opposite sides at frequencies with equal relative detuning 
from peak A and peak C, calculated using the TMM. We 
find strong nonreciprocity in all cases. The transmission 
in the two directions differs by up to a factor of 6 for 
peak A and a factor of 4 for peak C. Peak A also shows 
strong bistability for both directions. The switching in- 
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FIG. 5: PSTD results: Hysteresis curve of the structure for 
a Gaussian pulse incident from the left (a) and the right (b). 
The incident field has carrier frequency f c — 16.19 GHz (peak 
A), width parameter 1100// c , and peak amplitude 0.55. The 
arrows indicate the direction of change with time. 



tensity threshold for incidence from the left is almost 
an order of magnitude lower than for incidence from the 
right. With the given detuning from the linear transmis- 
sion peak frequency one finds no bistability for peak C 
but with stronger detuning bistability does appear. Such 
bistability in an optical diode was previously mentioned 
in 0,113. 

Finally, we simulate the transmission using the PSTD 
method. This time-domain method inherently includes 
the higher harmonics that may occur due to the nonlinear 
material. We use a spatial step-size of A x = Ao/75 and 
a temporal step-size of At — A x /(nc). The incident field 
was chosen as a long Gaussian pulse with peak amplitude 
0.55, carrier frequency f c — 16.19 GHz (peak A) and 
width parameter 1100// c , to operate in a regime where 
the envelope magnitude varies slowly. Figure shows 
the magnitude of the short-time Fourier-transform at the 
carrier frequency and at its third harmonic. Again, we 
find that the transmission strongly differs for incidence 
from the two sides. Transmission at f c is initially higher 
for incidence from the left but when the incident inten- 
sity reaches 0.04 transmission for incidence from the right 
becomes greater. On the other hand, the transmission at 
3/c is always much greater for the pulse incident from 
the left. At 3/ c the structure is essentially transparent, 
as the LH-material slabs have e r ~ fi r ~ n w 1 at that 
frequency. We find bistability for the carrier frequency, as 
well as, its third harmonic with incidence from either di- 
rection. The threshold switching intensities differ greatly 
for the two directions. The switching intensities for the 
two frequencies with a particular direction of incidence 
are very similar to one another. These results agree well 
with the TMM calculation. 

In our PSTD simulations we also encountered 
modulational-instability-like behavior in this system with 
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different threshold intensities for the onset of this insta- 
bility. This effect will be subject to further investigation. 

In conclusion, we have studied numerically the linear 
and nonlinear transmission of a novel type of multilayer 
structures composed of two different dielectric materi- 
als. We have shown that asymmetric periodic structures 



with nonlinear layers can demonstrate highly nonrecip- 
rocal transmission which is a major characteristic of an 
optical diode. We have analyzed the properties of the 
left-handed diode by employing the transfer-matrix ap- 
proach and direct pseudospectral time-domain numerical 
simulations, and have shown its bistable behavior. 
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